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Abstract 

We study fi — e conversion with sequential four generations. A large mass for the fourth generation 
neutrino can enhance the conversion rate by orders of magnitude. We compare constraints obtained from 
fi — e conversion using experimental bounds on various nuclei with those from fi — > e'y and \x — > eee. We 
find that the current bound from \i — e conversion with Au puts the most stringent constraint in this 
model. The relevant flavor changing parameter A^ e = V*^V e4 is constrained to be less than 1.6 x 10~ 5 for 
the fourth generation neutrino mass larger than 100 GeV. Implications for future fx — e conversion, [i — > e^f 
and [i —?• eee experiments are discussed. 
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I. INTRODUCTION 



The standard model (SM) of electroweak and strong interactions based on the gauge group 
SU(3)c x SU(2) L x U(l)y has been rigorously tested in many ways[l|, in different sectors involving 
gauge bosons, quarks and leptons. The flavor physics in the quark sector including CP violation 
is well described by the Cabibbo-Kobayashi-Maskawa (CKM) mixing matrixji), ^ in the charged 
current interactions of the W boson with three generations of quarks. Corresponding interactions in 
the leptonic sector are not on the same footing. In the simplest version of the SM, there are no right 
handed neutrinos, and the left handed neutrinos are massless. This theory conserves lepton flavor of 
each family separately, and no flavor changing neutral current (FCNC) processes are allowed in the 
leptonic sector. However, the observation of neutrino oscillations requires both a non-zero neutrino 
mass and lepton flavor violation (LFV). The SM has to be extended to give neutrino masses and 
explain their mixing through a mixing matrix in the charged currents analogous to the CKM matrix 
in the quark sector, the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) mixing matrix 0, There are 
many ways to extend the SM to have neutrino masses and mixing, for example, introduction of right- 
handed neutrinos to have Dirac mass, or to invoke seesaw mechanism to have Majorana mass, or to 
have loop induced masses. The existence of mixing among different lepton generations will induce 
FCNC interactions with LFV. To further understand the properties of the underlying theory, one 
should not only study whether a model can produce the correct neutrino masses and their mixing, 
it is also necessary to study possible implications for other FCNC interaction and confront them 
with existing and future experimental data. A particularly interesting process is \x — e conversion. 



/i — e conversion has been studied both theoretically and experimentally [6|-|9[. The quantity 
measuring the strength of this leptonic FCNC process is B A _^ e = ^convl^capt = + A(N, Z) — > 
e~ + A(N, Z))/r(/i~ + A(N, Z) — > v yi + A(N + 1, Z — 1). Experimental bounds on \i — e conversion 
for several nuclei have been obtained as B^ e < 7 x 10- 13 [Io|, B^ e < 7 x Hr u |ll|, B^ e < 
4.3 x 10- 12 [lij and B^ e < 4.6 x lO -11 j3] with the 90% c.f. level. These limits can give important 
information for leptonic FCNC interactions. 

There are many different ways to extend the SM to have neutrino masses and mixing, and 
therefore leptonic FCNC interactions. In this work, we study implications on \x — e conversion, in 
one of the simplest extension of the SM by introducing right-handed neutrinos to have Dirac masses. 
In this model, /i — e conversion cannot occur at tree level. At one loop level, it can be induced by the 
photon and Z penguin, and box diagrams. Since the Glashow-Iliopoulos-Maiani(GIM) mechanism 
is in effect, all leptonic FCNC processes involving charged leptons are proportional to ml./myy, here 
m u is the neutrino mass and m\y is the W boson mass. In the SM with three generations (SM3), 
due to the very tiny neutrino masses (< 0.2 eV), all LFV rates are extremely small. If the model is 
further extended to have a fourth generation (SM4), and if the fourth generation neutrino mass m U4 
is large enough, the factor m^./m^ can even become an enhancement factor, and therefore can have 
observable leptonic FCNC effects. In other words, the current experimental bounds can be used to 
constrain model parameters. FCNC interaction with SM4 has been studied for a long time 14| and 
can also help for solving some of open questions in FCNC quark interaction 15f. Leptonic FCNC 



effects, including fi — e conversion, in SM4 have also been studied extensively |l6Ml9|. In our work 



we carry out a more detailed systematic study by taking account of all available \i — e conversion 
experimental results. We compare constraints obtained from \x — e conversion using experimental 
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bounds on various nuclei with those from /i — > cy and /i — > eee. We find that the current bound 
from fi — e conversion with Au puts the most stringent constraints on model parameters. The 
relevant flavor changing parameter A^ e = V*^/ e i is constrained to be less than 1.6 x 10~ 5 for m U4 
larger than 100 GeV. Future improved experiments on [i — e conversion, /i — > c-f and /i — > eee will 
further constrain model parameters. 



II. THE FCNC INTERACTIONS FOR fx-e CONVERSION IN SM4 

The minimal set of particle contents in the SM is the gauge bosons of the gauge group 
SU(3) C x SU(2) L x U(1) Y with couplings g s , g and g', the Higgs doublet $ : (1,2,+1) with the 
component fields and vacuum expectation value (vev) v = (h) given by $ = (co + , (v + h + iz)/\/2) T ', 
the three generations of left and right handed up and down quarks, Ql '■ (3, 2, 1/3), Ur : (3, 1,4/3), 
D R : (3, 1, —2/3), and the three generations of left handed leptons L L : (1, 2, —1) and right handed 
charged leptons, Er : (1, 1, —2). We use the hypercharge normalization that Q = I$ + Y/2. Depend- 
ing on the nature of the neutrinos, there may be right handed neutrinos, Nr : (1,1,0) which may 
pair up with the left handed neutrinos to have Dirac or by themselves to have Majorana masses. 

For a model with Dirac neutrino mass only, the charged lepton and neutrino masses are generated 
by the Yukawa interactions 

L = -L L Y E ^E R - L L Y N $N R + h.c. (1) 

where $ = ia 2 ^*, E R = (e, /i, r, ...)^, N R = (z/ e , z/ T , ...)^ for n generations of leptons. Y E and Y v 
are n x n matrices. 

After the Higgs boson develops a non-zero vev, the leptons will obtain their masses with the 
mass matrices for charged leptons and neutrinos which are given by 

M E = Y E ^=, M U = Y N ^=. (2) 

In the basis where the mass matrices are all diagonal, the flavor changing interactions involving 
leptons are contained in the PMNS mixing matrix U = (Vij) defined by 

L = --^E L rUv L W; + h.c. (3) 

where El = (e, /i, r, ■■■)l, vl = [yi, ^2, ^3, an ^ U \s a. n x n unitary matrix for n generations. 
For SM3, U can be parameterized as[l| 

(C12C13 S12C13 s^e lS \ 

-S12C23 - Ci2S23Sl 3 e t ' 5 C12C23 - S 12 S 2 3S 13 e' 15 S23C13 , (4) 
S12S23 - Ci2C23Si 3 e J<5 -C12S23 - S 12 C 2 3Sl3e lS C23C13 / 

where Cy = cos^- and = sin^- with 8^ being rotation angles and 5 being the CP violating 
phase. 

The charged lepton masses have been measured with great precision. But the neutrino masses 
are not. The absolute values for neutrino mass scale, although not known, are constrained to be 



3 



small (m„ < 2.3 eV by tritium beta decay experiment 20j, and Yl m i < 0-61 g V by cosmological 
observations (2 11] ) . and the mass differences are known from neutrino oscillation data. The mixing 
angles are also constrained by various experimental data in particular from neutrino oscillation 
data. We list these constraints on neutrino mixing and masses in the following 



sin 2 (2# 12 ) = 0.87 ±0.23 (22|, sin 2 (2# 2 3) > 0.92 |23|, sin 2 (2#i 3 ) < 0.15 CL 

(2.43 ±0.13) x 10- 3 eV 2 |25 



90% 24 



Am 2 21 = (7.59 ± 0.20) x 10~ 5 



22] 



\ Am li\ 



(5) 



where Am?- 



mf — m 2 , and rrii(i = 1-3) are the mass eigenvalues of neutrinos. Recently T2K 
has obtained new evidence j26|] showing that the angle #13 may indeed be non-zero with the 90% c.l. 
allowed region 0.03(0.04) ~ 0.28(0.34) for sin 2 (2# 13 ) for normal (inverted) neutrino mass hierarchy. 

In this model FCNC processes do not occur at the tree level, but it can happen at one loop 
level. Due to GIM mechanism of the model for lepton sector, amplitudes for any leptonic FCNC 
processes with charged leptons are proportional to Am-j/m^. Since squared differences of neutrino 
masses are constrained as in above, all FCNC processes involve charged leptons including /i — e 
conversion, are very small which are well below experimental sensitivity. To have observable FCNC 
effects, further extensions are needed. A simple way which can overcome the problem of small 
Amfj/rriw ratio is to have a heavy neutrino in the model. This can be achieved for a model with 
four generations (SM4). 

The existence of a fourth generation of quarks and leptons is not ruled out although there are 
severe constraints. The heavy neutrino must be heavier than m z /2 such that the precisely mea- 
sured Z decay width would not be upset. The LEP II data constraint the mass to be even larger 

The constraint from electroweak precision data can also be 

Perturbative unitar- 
Critical 



with m Vi > 90.3GeV at 95%CL (27 



evaded by imposing appropriate mass splitting between fourth generations 28 



ity bound on a heavy fourth generation Yukawa coupling constant has been considered [29 



masses determined by perturbativity for charged leptons and neutrinos are about 1.2 TeV 29[. The 



fourth generation neutrinos, if exists, with a large mass in the range about 100 GeV to the unitarity 
bound are not ruled out. With the addition of the fourth generation, the PMNS matrix U will 
become a 4 x 4 unitary matrix. If mixing matrix elements in Vn are not extremely small, larger 
FCNC effects can be generated at one loop level. From the neutrino oscillation data, the constraints 
on the mixing of the fourth generation with the first three light neutrinos are indeed allowed to be 
not extremely small as the large error bars in eq.(j5]) indicate. There is a chance for large FCNC 
effects in the lepton sector. 

The strength of muonic FCNC processes is governed by powers of A Me = V* A V eA . Loop induced 
FCNC interactions relevant to /1 — e conversion can be obtained by standard loop calculations. We 
present the one loop effective Lagrangian in two terms in the following, in the limit of neglecting 
terms suppressed by light neutrino mass squared divided by the W mass squared, 



L{n erf) 



L{\± ->• e + qq) 



-^=-^-^Xf l eG 2 (x 4 )ea' J ' u (m fl R + m e L)/j,F^ , 
Gf e 2 

= -—j=— [V u (x 4 )u'y II u + V d (xf)djfj,d 
±A u (x 4 )«7 M 75W + A d (x A )d-i^ b d] ej^Lfi , 



(6) 
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where o> = i[7 At ,7„]/2, i?(L) = (1 ± 7s)/2, x 4 = m 2 v Jm 2 Wl and 



G x {x) 

G 2 (x) 

G z (x) 
G B (x) 



-CPGi{x) + -^(/| - Q q s 2 w )G z (x) 
Zs w 



— -hGz(x) + 
s w 

x{\2 + x - 7x 2 ) 

12(x - l) 3 
x(l — 5x — 2x 2 ) 



2 ' 

_ x 2 (12 - lOx + x 2 
6(x - l) 4 

3x 3 



lnx 



4(x - l) 3 2(x - l) 4 

x(x 2 - 7x + 6) x(2 + 3x) 



4(x - l) 2 



x 



x 



[X 



1) : 



4(x - l) 2 
lnx, G B (x 



In a; , 
lnx , 



4G B (x) 



GUx) 



Gr(x) . 



(7) 



In the above, Q q is the electric charge of the quark q. 

We also present the effective Lagrangian for /i — > eee in the following which will be used when 
we compare constraints on parameters, 

G F e 2 



L{n — > eee) 



y/2 4vr 2 



A 



fj,e 



G 2 (x 4 ) erj^e eia^{m^R + m e L)\i 



+ (a L (x 4 )L + a R (x A )R) e e^^Lfi 



where 



a L (x) = GAx) + 



1 

2 



Gz(x) - 



2s 2 



G B {x) , a R (x) = Gi(x) + G z {x) 



(8) 



(9) 



The term L(fi — > ej) is generated by the electromagnetic penguin. The loop generates the 
function G 2 (x) for dipole operator. For the terms L(fi — > e + gg) and L(/i — > eee), there are several 
contributions. Connecting the photon in the dipole operator, this generates the terms proportional 
to G 2 (x). The electromagnetic penguin also generates a charge radius term. With the photon 
connected to a quark or an electron, the terms proportional to G\{x) are generated. Replacing the 
photon by a Z boson, terms proportional to Gz(x) are generated. Finally, the box diagram produce 
the terms proportional to Gb{x). We note that in the large x limit the function Gz{x) linearly 



increases, while others are not. The results above agree with those in Ref. 18 



III. n ->■ e CONVERSION CONSTRAINTS ON MODEL PARAMETERS 

To obtain \i — y e conversion rates from the FCNC interaction discussed in the previous section, 
one needs to convert quarks into hadrons which depend on hadronic matrix element calculations. 
Theoretical efforts have been made by several group to calculated these matrix element sjo)-^. We 
will use the results in Ref. for a consistent calculation. In this framework there are several quark 



5 



level contributions which are parameterized as the following effective Lagrangian 



[m^ u (A R R + A L L)fiF^ + h.c] 
G 

-j= [e(9LS(g)R + 9RS( q )L)fi qq + e(g LP[q) R + g RP{q) L)ii q^q + h.c] 
G 

-j= [e(gLV( q )^L + g R v( q )l^R)^qi^q + e(g LA { q )l^L + g RA{q) ^R)iiqy^ 5 q + h.c] 



Gf 
V2 



7: e (9LT{g)^R + gRTiq^ty/iqa^q + h.c. 



(10) 



In SM4, at the one loop level, we have, g(L,R)(S,p,T)( q ) = 0, gnv( q ) = 0, and 



2 

Ar = g^^ 2 A Me C 2 (x4) , A L = —Ar , 9LV{q) = -^^iieVq ■ (H) 

The contribution from Al can be neglected compared with that from A R . The contribution from 
9l(r)a terms proportional to A q in eq.((6]) is suppressed because the fraction of the coherent process 
is generally larger than that of the incoherent one approximately by a factor of the mass number 
of the target nuclei jsj]. 

The quantity B A e measuring the leptonic FCNC effect in ji — > e conversion to the leading order 
in SM4 is proportional to 

\A R D+~g^V^ + ~gtiv^\\ (12) 

where 

9lv = 2 9lv(u) + 9LV{d) , 9lv = 9lv(u) + 2 9LV(d) • (13) 
Combining these, we obtain 

B A 



ArD(A) A r D(A) 



(14) 



where 



R° U .JA) = ° f T\ \D(A)\ 2 . (15) 

^ eK 1 i927r 2 r4 P i v n y 1 

The relevant parameters are listed in Table [J which are evaluated using method I in Ref.jsj]. 

In FigJT] we show constraints on A Me as a function of the fourth generation neutrino mass m U4 
using experimental upper limits on \x — e conversion for Au, S, Ti and Pb nuclei. The most stringent 
constraint comes from Au. We see that A Me is constrained to be less than 1.6 x 10~ 5 for m Vi larger 
than 100 GeV. Such a small value for is well within the error bars allowed ranges from neutrino 
oscillation data. To see whether fx — e conversion provides the strongest constrains on the model 
parameters, one needs to check with other leptonic FCNC processes. In the next section we compare 
the fx — e conversion constraints with those obtained from other two well known leptonic FCNC 
processes fx — > ey and fx — > eee. 
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A 


D(A) 


yW(A) 


yW)(A) 




27 A1 

32 c 
16 D 

iTi 

79 7 AU 

Ifpb 


0.0524 
0.0864 
0.189 
0.161 


n ni fii 
0.0236 
0.0396 
0.0974 
0.0834 


0.0236 
0.0468 
0.146 
0.128 


0.0028 
0.0041 
0.0039 
0.0027 



TABLE I: The relevant parameters for /x — e conversion processes, which are evaluated by using method I 
in Ref.fl. 
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FIG. 1: Constraints on |Au e | as a function of neutrino mass m Vi from current experimental bounds on fi— e 
conversion, /x — > cy and /x — > eee. 



IV. COMPARISONS WITH CONSTRAINTS FROM xx -> ej AND /x -> eee 



xx — >■ e7 and /x — > eee are two well known processes which provide stringent constraints on lep- 
tonic FCNC interactions. Experimental measurement for /x — >■ cy is not an easy task. Nevertheless, 
several groups have made efforts in obtaining limit on this decay. No positive result on the ob- 
servation of this decay has been reported. The best upper bound at the 90%CL, set more than 

10 years ago, is b(? -> e 7) < i.2 x io-»0. Recen ; ly MEG collaboration has reported their 

preliminary result on the upper limit with B(fi — > ej) < 1.5 x 10 -11 based on their 2009 physics 



data collection |31|. This bound is slightly weaker than the best upper bound. It is, however, in- 
triguing to note that they did have a few events well separated from backgrounds in their photon 
and electron energy cuts plot. However, due to low statistics, no positive claims can be made. 



MEG experiment has potential to improve the sensitivity by two orders of mag nitudes, 10- 13 |32j 
Experimental studies of /x — > eee has put a stringent limit on the branching ratio. The current best 
limit is B(fi -»■ eee) < 1.0 x 1O" 12 0- 

Using the effective Lagrangian L(/x — > &y) in eq. ([6]), one can obtain the branching ratio B(fi — > 
ej) = — > ej)/r(fi — > evv). Neglecting m e , one obtains 

fi(/x^e 7 ) = 384tt 2 (\A l \ 2 + \A r \ 2 ) , (16) 

where Ar l are defined in eq.liTTl). If the mass m e is kept, one should divide, in the above expression, 
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FIG. 2: Ratios of B£_ e \/B(fi — > &y) and B(fi — > eee)/B(fi — > ej) as functions of m h 



a phase factor I{x) 



ix- 



sx 



-X 



■12x 2 lax with x = m 2 e /rn? . Radiative corrections from QED 



also modify the above expression by dividing a factor 1 + Sqed with 5qed = (a/27r)(25/4 — 7r 2 
The branching ratio for \i — > eee can be deduced from L(/i — > eee) in eq.(jHJ). We have 



B(fi — > eee) 



167T 2 ' ^ 



[a 2 R (x 4 ) + 2a 2 L (xi) 



- 4G 2 (x4)(a fi (x 4 ) + 2a L (x 4 )) + 4G^(x 4 ) (4 In - — 



(17) 



The constraints on the parameter A Me as functions of m„ 4 for bounds from fi — y ej and \x — > eee 
are also shown in Fig{T] The constraint obtained from \i — > eee is weaker for neutrino mass m u4: 
less than 255 GeV than that from /i — > ej, but becomes stronger for larger mass for the current 
experimental bounds. This property is mainly caused by Z-penguin contributions. However, the 
constraints from \x — e conversion on Ti and Au are better than those from fi — > e^ and \i — > eee. 
The best constraint comes from \i — e conversion on Au nuclei. 

In FigJ2]we show the ratios of B^_ e JB{^ — > ej) and B(ji — > eee)/B(n —> ej) as functions of 
m U4 . For a given fourth generation neutrino mass, the ratios are fixed independent the value of 
A Me . We see that in the region with m VA larger than 100 GeV, the ratios for \x — e conversion of 
various nuclei are larger than one implying that if the experimental values for the fi — e conversion 
rate and \i — > e^ branching ratio are similar, the constraint from \i — e conversion will give stronger 
constraints on the parameter A Me which explains the fact that the current experimental bounds on 
fi — e conversion from Ti and Au provide stronger constraints. When m„ 4 becomes larger ratios 
become larger due to the fact that the Z penguin whose contribution increases with m^ 4 for \i — e 
conversion and [i — > eee, but fi — > e^ increases only logarithmically in large m„ 4 limit. 

If one takes the well separated events of the MEG data as \i — > ej events, the branching ratio 
implied would be of order 10 -12 , we will take an assumed branching ratio of 3 x 10 -12 to show the 
implication on the parameters required to produce it. We show the results in FigJSl We see that such 
a possibility is ruled out by fi — e conversion on Au nuclei in this model (see also Figfl]). In FigJHl 
we also plot possible constraints on parameters using several projected experimental sensitivities 
on the fi — e conversion, fi — )■ e^ and fi — > eee. We see that the improved experiment MEG on 
fi — > e7 and MuSIC[34| on fi — > eee can obtain better constraints than the current Au bound. 
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FIG. 3: Constraints on 4th-generation mixing coupling |A« e | as a function of neutrino mass m Ui from 
several future experimental sensitivities on fi — e conversion, ft — )■ and fi — > eee. 



Future experiments Mu2E 35J /COMET [36j and PRISM 37J for fi — e conversion using Al and Ti, 
respectively, can obtain much better constraints. 



V. CONCLUSIONS 



We have studied fi—e conversion in the SM with four generations. A large mass m U4 for the fourth 
generation neutrino can enhance the conversion rate by orders of magnitude. We have compared 
constraints obtained from \i — e conversion using experimental bounds on various nuclei with those 
from fx — > e7 and fx — > eee. We found that the current bound from fx — e conversion with Au puts 
the most stringent constraint in this model. The relevant flavor changing parameter A Me = V^* 4 VJ=4 is 
constrained to be less than 1.6 x 10~ 5 for m Ui larger than 100 GeV. The goal of MEG experiment is 
to have a sensitivity of 10~ 13 for B(/j, — > ej), which can compete that of fi — e conversion in Au in 
this model. Therefore, the model prediction can be tested bv combining these coming data. There 



are several proposed experiment for fi — e conversion, Mu2E 35[, COMET 36j and PRISM 37|, which 



will tell us detailed information of lepton flavor mixing. In addition, successful running and energy 
upgrade of Large Hadron Collider will search for sequential fourth generation fermions. The fourth 
generation model will be explored directly and indirectly by future experiments in great detail. 
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